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^ ! Abstract 

^ I We prove an optimal 4G Theorem for the Gaussian kernel. We also 

propose a new general method of estimating Schrodinger perturbations 
^ of transition densities, and give applications to the Gaussian kernel. 

1 Introduction 

! A Schrodinger perturbation is an addition of an operator of multiplication 

r~| ' to a given operator. On the level of inverse operators, the addition results 

■ in resolvent or Duhamel's or perturbation formula, and under certain con- 

ditions it yields von Neumann or perturbation series for the inverse of the 
perturbation. The subject is very wide, and we intend to study it in the case 
when the inverse operator is an evolution semigroup, in fact, a transition 
J> . density. In this case a convenient and simple setup is that of integral opera- 

I tors on space-time, and the perturbation series has an exponential flavor due 

\^ ' to repeated integrations on time simplexes. We propose a general method for 

'sj" . pointwise estimates of the series, and we show its versatility by estimating 

transition densities of Schrodinger perturbations of heat kernels on R'^. 
' In an earlier work, Bogdan, Jakubowski and Sydor [6] developed a tech- 

nique for sharp pointwise estimates of Schrodinger perturbations p of tran- 
sition densities p and more general integral kernels on state space X by 
functions q > 0. The method rests on the assumption 



p{s, X, u, z)q{u, z)p{u, z, t, y)dzdu <[ri + Q{s, t)]p{s, x, t, y), 



X 
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where s < t, x,y £ X, rj < oo, and < Q{s,u) + Q{u,t) < Q{s,t). The 
left-hand side of (1) defines the term pi{s,x,t,y) in the perturbation series, 



oo 

p{s,x,t,y) = '^pn{s,x,t,y) 

n=0 



(2) 



(see Section 2), and so pi{s, x,t,y)/p{s, x,t,y) < r] + Q{s,t). The bound is 
uniform in space and locally uniform in time, and it propagates as follows. 



Pn{s,x,t,y) < pn-i{s,x,t,y) 

n 

< p{s,x,t,y) Y[ 



k=l 



v + 



Qis,t) 

n 

Q{s,t) 

k 



(3) 



Furthermore, if < ?] < 1, then (2) and (3) yield. 



p{s,x,t,y) < p{s,x,t,y) 



1 — 7] 



and if = 0, then 



p{s,x,t,y) <p{s,x,t,y)e^^'^^l 



(4) 



(5) 



The above estimates are sharp i.e. the ratio of upper bound of p and (the 
trivial lower bound) p is bounded locally in time. In fact, as shown by 
Bogdan, Hansen and Jakubowski [4, Example 4.3 and 4.5] the exponen- 
tial factors in (4) or (5) are nearly optimal. The estimates also apply to 
very general integral kernels on space-time, without assuming Chapman- 
Kolmogorov equations. It is now crucial to verify (1) for given p and q. To 
this end we usually try to split and estimate the singularities of the inte- 
grand {u,z) I— 7- p{s,x,u,z)q{u,z)p{u, z,t,y) in (1). This is straightforward 
if q satisfies a suitable Kato-type condition and 3G Theorem holds for p 
(see [6, Remark 2] and the discussion in Section 2 and Section 4 below). 
The latter is the case, e.g., for the transition density of the fractional Lapla- 
cian as described by Bogdan, Hansen and Jakubowski [3, Corollary 11], [4, 
Example 4.13], and for the potential kernel of the stable subordinator [6, 
Example 2], because the functions have power-type asymptotics. 

However, 3G fails for the Gaussian transition density due to its exponen- 
tial decay. Motivated by this obstacle, the results of Zhang [26, 27] and the 
arguments of Jakubowski and Szczypkowski [17, 18], in Section 2 we propose 
a more flexible method of estimating Schrodinger perturbations of transition 
densities. The method employs an auxiliary transition density p* , as an 
approximate majorant of p substituting for p{-,-,t,y) in (1). In the second 
part of the paper we test our method on Gaussian-type estimates. Namely, 
in Section 3 we elaborate [26, inequality (4.4)] by giving a sharp estimate 
involving four Gaussian kernels (hence 4G), one of which may be interpreted 
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as the majorant. Then, in Section 4 we obtain precise (but necessarily not 
sharp) Gaussian estimates for the fundamental solution of Schrodinger per- 
turbations of second order parabolic differential operators, recovering and 
improving existing results, which we discuss there at some length. We also 
describe a connection to second order differential operators. In Section 5 we 
give miscellaneous methodological comments. 

Our inspiration mainly comes from [18] and [6]. Our ideas are also similar 
to those developed for Gaussian estimates in [26]. In particular, the condi- 
tion (16) for our main Theorem 1 may be considered as a generalization of 
the Main Lemma 4.1 of [26], while the 4G inequality in Theorem 2 yields 
an alternative, synthetic justification of that lemma. Furthermore, the proof 
of [26, inequality (4.4)], yields (24), except for the optimal constant M. It 
is thus of interest that the approach of [26], which was taylor-made for the 
Gaussian kernel, has a more general context given by Theorem 1. Note- 
worthy, our approach is not restricted to Gaussian-type kernels. Further 
applications, e.g. to perturbations of the transition density of the 1/2-stable 
subordinator, will be given in a forthcoming paper. 



2 General estimates of transition densities 

Let X be an arbitrary set with a cx-algebra M and a (non-negative) cr-finite 
measure m defined on M. To simplify the notation we will write dz for 
m{dz) in what follows. We also consider the Borel subsets B of M, and the 
Lebesgue measure, du, defined on M. The space-time, MxX, will be equipped 
with the cr-algebra B x M and the product measure dudz = dum{dz). 

We will consider a measurable transition density p on space-time, i.e. we 
assume that ptMxXxMxX— > [0,oo] is B x M x B x 7W-measurable and 
the Chapman-Kolmogorov equations hold for all x,y £ X and s < u < t: 

J p{s, X, u, z)p{u, z, t, y) dz = p{s, x, t, y) . (6) 
X 

Let g : M X X — >• [0, oo] be (nonnegative and) B x A^-measurable. The 
Schrodinger perturbation p of p by g is defined by the series (2), where 
Po{s,x,t,y) =p{s,x,t,y), 

t 

pi{s, x,t,y) = j j p{s, X, u, z)q{u, z)p{u, z, t, y)dzdu, 

s X 
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and for n = 2,3, . . ., 

t t t 

Pn{s,x,t,y)= /••• / / pis,x,ui,zi)q{ui,zi) 



S Ul U„-l 

p{ui, zi,U2, 2:2) •• • q{un, Zn)p{un, Zn, t, y)dzn ■ ■ ■ dzidu^ ■ ■ ■ dui- (7) 
By Fubini-Tonelli, for n = 1, 2, . . ., we have 

i 

Pn{s,x,t,y) = J j p{s,x,u,z)q{u,z)pn-i{u,z,t,y)dzdu , (8) 

s X 

and 

t 

Pn{s,x,t,y) = j j pn-i{s,x,u,z)q{u,z)p{u,z,t,y)dzdu. (9) 

s X 

By [3, Lemma 1], for all s < n < t, x, y G X and n € Nq = {0, 1, . . .}, 

n „ 

Pm{s,X,U,z)pn-m{u,Z,t,y)dz = Pn{s,X,t,y) . (10) 

By [3, Lemma 2], Chapman- Kolmogorov equations hold for p. Clearly, p > p. 

Remark 1. The perturbation, say pv, is given by the same formulae if 
y : R X M"^ — >• C (takes on complex values), provided is finite. Then pv 
converges absolutely and 

\pv\ <P\v\- (11) 

A detailed discussion of signed real-valued perturbations of transition den- 
sities is given in [3], with a positive lower bound for pv which follows from 
Jensen's inequality. A probabilistic interpretation of pn and pv may also be 
found in [3]. 

Below we focus on upper bounds oi p = pq for transition densities p and 
g > 0. This is a less general setting than that of [6], but within this setting 
our bounds, (17) and (18) below, hold under more flexible condition (16) on 
p and q. Namely, we consider another measurable transition density p* and 
C > 1 such that for all x,y ^ X and s < t, 

p{s,x,t,y) <Cp*{s,x,t,y) . (12) 

We can estimate the cumulative effect of the integrations involved in (8). 
The following result is an analogue of [18, Lemma 5] and [4, Example 4.5]. 
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Lemma 2. Let 9 > and sq < . . . < Sk = t be such that 

Si + l 

j j p{s,x,u,z)q{u,z)p*{u,z,Si+i,y)dzdu <ep*{s,x,Si^i,y), (13) 

s X 

for all i = 0, . . . ,k — 1, s £ [si, Si+i] and x,y £ X . Then for every n G Nq, 



Pn{s,x,t,y) < y yC^p*{s,x,t,y), s £ [sq, si\, x,y £ X. (14) 

Proof. For k = 1, the estimate holds for n = by (12), and then it holds for 
all n > 1 by induction, (8) and (13): 



Pn{s,x,t,y) = J J p{s,x,u,z)q{u,z)pn-i{u,z,t,y)dzdu 

s X 

t 

an— 1 



< 0" ^ j j PI'S] 2;, n, z)(7(n, z)p*(u, t, y)d2;(iu 

s X 

< 9"'Cp*{s,x,t,y), where sG[so;Si], x,y £ X. 
li k > 2, then by (10), induction and Chapman- Kolmogorov for p* , 

n „ 

Pn{s,X,t,y) = ^ / Pmis,X,Sk-l,z)pn-niisk-l,Z,t,y)dz 

/ ("'^'^~'^)o"'C''~Vis,x,Sk-i,z)e"-"'Cp*{sk^i,z,t,y)dz 



n + k 1 ^-"(^^(^^^.^i^y)^ [is£[so,si],x,yGX,n£'No 



k-1 

□ 

In passing we note that the assumption and conckision in the statement 
of [18, Lemma 5] need a slight strengthening for the induction to work prop- 
erly: each tj+i in the assumption there should be replaced by r in 
and each t in the conclusion should be replaced by r in (then one 

proceeds as in the proof of Lemma 2 above) . The correction does not influ- 
ence applications of Lemma 5 or other results in [18]. 

We further let Q: M x M — )• [0, oo) be regular superadditive, meaning that 

Q{s,u) + Q{u,t) < Q{s,t) , if s<u<t, (15) 

Q{s,t) = if s > s I— )• Q{s,t) is right-continuous and t i— )• Q{s,t) is left- 
continuous (the continuity assumptions are rather innocuous, as we explain 
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later on in Lemma 12 and Lemma 13). We see that t i— t- Q{s,t) is non- 
decreasing and s I— 7- Q{s,t) is non-increasing. For instance, if /i is a Radon 
measure on M, then Q(s,t) = € M : s < u < t}) is regular superaddi- 

tive. A regular superadditive Q is infinitely decomposable in the following 
sense. 

Lemma 3. Let s < t, k ^ N and 6 > be such that Q{s,t) < kO. Then 
s = sq < si < . . . < Sk = t exist such that Q{si^i, Sj) < 9 for i = 1, . . . ,k. 

Proof. We may and do assume that k > 1 and (k — 1)6 < Q{s,t) < kO. 
Let Si = m.i{u : Q{s,u) > i9} for i = 1, . . . ,k — 1. If s < u < Si, then 
Q{s,u) < i9, and so Q{s,Si) < iO. li si < u < Sj+i, then Q{s,u) > i9 and 
Q{s,u) +Q{u,Si^i) < Q{s,Si+i) < {i + 1)9, thus Q{u,Si+i) < 9. Letting 
u ^ Si we obtain Q{si, Sj+i) < 9, which is also true if Si = Sj+i. □ 

Definition 4. We write q £ J\f{p,p* ,C,rj,Q) if q > is measurable on 
space-time, p and p* are measurable transition densities, C >1, r]>0, Q is 
regular superadditive, and (12) and the following inequality hold for all s < t 
and x,y £ X , 



p{s,x,u,z)q{u,z)p*{u,z,t,y)dzdu< [q + Q{s,t)]p*{s,x,t,y) . (16) 



s X 



The terms r] and Q{s,t) of (16) propagate differently in estimates of p„ 
below. We may think about rj as giving a bound for instantaneous growth 
of mass, while Q gives a cap for growth accumulated over time (see [6] and 
[4] for such insight). 

Theorem 1. If q £ M{p,p*, C, rj, Q), then for all s < t and x,y £ X, 



p{s,x,t,y)<p*is,x,t,y)\^-:^—^j , if < < 1/2, (17) 

and 

p{s,x,t,y)<p*{s,x,t,y){2Cf+^'^^'''\ if r/ = 0. (18) 

Proof. Let A; G N. By Lemma ?>, sq = s < si < . . . < Sk = t exist such that 
Q{si-i, Si) < Q{s, t) /k for i = 1, . . . ,k. By Lemma 2 with 9 = r]-\- Q{s, t) /k, 
and by Taylor's expansion, for all x,y £ X we get 

oo 

p{s,x,t,y) < ^ \pn{s,x,t,y)\ 

n=0 



n=0 

k 



^ ' P*{s,x,t,y), 



l-T]-Q{s,t)/k 
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For > 0, we choose A; € N such that {k — l)rj < Q{s,t) < krj. For rj = 0, 
we choose k satisfying {k — l)/2 < Q{s,t) < k/2. This ends the proof. □ 

Analogous results hold if we state our assumptions and conclusions for 
in a finite time horizon: — oo <ti<s<t<t2< oo. If Q in Theo- 
rem 1 is bounded, then p < const, p* uniformly in time. We may consider 
q G M(j),p* ,C,r],Q) with r] < 1/2, qi{u,z) = g(n, z) l[o.i](n), and bounded 
superadditive function Qi{s, t) = Q{{s V 0) A 1, (t A 1) V 0). Then, 



p{s, x,u, z)qi{u, z)p* {u, z,t,y) dzdu <[r] + Qi{s,t)]p* {s, x,t,y) . 



s X 



Thus, by Theorem 1, p < cp* uniformly in time. If p* is not comparable 
with p in space then the estimates in Theorem 1 cannot be sharp. This is 
regrettable, but quite common, e.g., in Schrodinger perturbations of Gaus- 
sian kernel discussed later on in the paper. The role of p* is similar to that 
of / in [4, Theorem 3.2], but the results of [4] do not apply in the present 
setting, if p ^ p* . If (16) holds with p* replaced by p, then we may take 
p* = p and C = 1 in (12) and Theorem 1. However, in this case a more 
efficient inductive argument of [6] gives better estimates (4) and (5) above. 
If q{u, z) < f{u), then we may take Q{s, t) = f{u)du, r] = and p = p* . 
In fact, if q{u,z) = f{u), then pn{s,x,t,y) = Q{s,t)'^p{s,x,t,y)/n\ and 
p{s, X, t, y) = eQ('''*)p(s, X, t, y). 

Less trivial applications of Theorem 1 require detailed study of p and 
g, and judicious choice of p* . In particular, to estimate p for given p, p* 
and q, we wish to verify (16). This task may be facilitated by splitting the 
singularities of p and p* in the integral of (16). Various versions of the 3G 
Theorem are used to this end, see [3], Bogdan and Jakubowski [5], and the 
"elliptic case" of Cranston, Fabes and Zhao [9], and Hansen [14]. For instance 
the transition density of the fractional Laplacian enjoys the following 3G 

p{s,x,u,z) Ap{u,z,t,y) < cp{s,x,t,y), 

where x,z,yGM'^, s < u <t ([5, Theorem 4]), and this yields 

p{s, X, u, z)p{u, z, t, y) = [p{s, X, n, z) V p{u, z, t, y)] [p{s, x, u, z) A p{u, z, t, y)] 
< c [p{s, X, n, z) + p{u, z, t, y)]p{s, x, t, y). 

In this situation we can use p* = p m. Theorem 1 to estimate jj, provided 
t t 

p{s,x,u, z)q{u, z) dzdu + c / / q{u, z)p{u, z,t,y) dzdu < r] + Q{s,t) . 



s X s X 



Noteworthy, 3G fails for Gaussian kernels, and for such kernels the methods 
of [6] fall short of optimal known results. This circumstance largely motivates 
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the present development. Below we show how to estimate quite general 
Schrodinger perturbations of Gaussian kernels by means of Theorem 1 . The 
application depends on a 4G estimate, which is stated in (24) of Theorem 2 
and partially substitutes for 3G. We note that Theorem 2 improves [26, 
(4.4)], since we give an optimal constant in (24). Explicit constants matter 
in applications, because we specifically require rj < 1/2 in Theorem 1. 

3 Estimates of the Gaussian kernel 

As usual a\/ b = max{a, b} and a f\b = min{a, 6}. Let < a < oo, and 

(19) 



L{a) = max 

T>aVl/a 



max 

T>aVl/a 



T — OL 

ln(l +r) In(ar) 

1 + T 



In ( 1 H — ) — In a + — — In(aT) 
T 1 + T 



Clearly, L{a) < oo, and r = a V 1/a yields L{a) > ln(l + a V 1/a). By an 
application of calculus, L{a) = ln(l + a) if (and only if) a > e^l'^ . We let 

/(r,x) = In r + x^/r, r > 0, x > 0. 

Lemma 5. If a > 0, L = L{a), ^,7? > 0, and t > 0, then 

/(I + r, e + r/) < /(1, V /(ar, r?) + ^ + L. (20) 

r 

If L < L[a), then the inequality fails for some £,,r] >0 and r > 0. 
Proof. We first prove the following implication: 

If ^+ln(ar)<^^ then ln(l + r) < -^^^^^^ + L . (21) 
ar r(l + r) 

To this end we consider two special cases: 

Case 1. rj^ /{ar) + In(ar) < and r] = t^^ , 

Case 2. rj^ /{ar) + In(ar) = and rj < t^, . 

Case 1 implies that (r/a — 1).^^ + ln(Q;r) < 0. This is possible only if 
r < a V 1/a, whence ln(l + r) < ln(l + a V 1/a) < L{a), which verifies (21). 

In Case 2, if r < a V 1/a, then ln(l + r) < ln(l + a V 1/a) < L{a) again. 
For r>aVl/awe consider C = ?(^) as a function of 77, and we have 

c^ir^) := ln(l + r) - < L. 

r(l + r) 

Indeed, we see that the condition rj < holds automatically since rf /t"^ < 
rf' /{ar) = ^"^ — In(ar) < Our assumption now reads = rf /{aT) + 
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In(aT), where (,,rj > and r > a V 1/a. Thus, ^' = r}/[aT^). Note that 
0(0) = ln(l + r) — r ln(ar)/(l + r) < L. Furthermore, 

</>'(r/) = -2(l + r)-i(re-r?)(e'-l/r). 

We have (j)'{r]) = only if ^' = 1/r, or = rj/a, and then r/^/ (ar) +ln(ar) = 
rf /a^ and i;^>(r/) = ln(l + r) — (r — a) ln(ar)/(l + r). This in fact shows that 
L = L{a) is sharp in (21), see (19). Furthermore, (p [rj) < if ^' > 1/r, or 
(r/a — l)rf /a > rln(ar), in particular if rj is large. Thus, (/) is decreasing 
for large r], which yields (21) in Case 2. 

Consider general ^,r] and r > in (21). If r/ > r^, then decreasing r] to 
strengthens (21), so eventually we are done by Case 1. If 77 < r.^, then we 
increase r] and strengthen the consequent of (21), getting under Case 1 or 2. 

Putting (21) differently, ln(l + r) + (^ + 7?)2/(l + r) < i^ + rf/r + L, pro- 
vided "rf /{ar) + In(ar) < Therefore we have (20) under the assumption 
f{aT,r}) < /(l,^), and the constant L cannot be improved. In particu- 
lar, (20) holds if f = f{aT,r}). Decreasing ^ keeps (20) valid because 
/(I + r, + 77) then decreases, too. □ 

For d G N, a > 0, s < t, and x, y G M'^, we consider the Gaussian kernel 

ga{s,x,t,y) = [A7T{t-s)/a]-'^/^e^p{-\y-x\y[4{t-s)/a]}. (22) 

Let g = gi- Clearly, if < a < 6 < 00, then for all s < t and x, y G M*^, 

gb{s, X, t, y) < {h/af'^ga{s, x, t, y) . (23) 

We note that 3G fails for ga, because 

ga{G,Q,t,y)^ga{t,y,2t,2y) ^ (47rVa)-'^/^e-l^l'/("^/°) ^ ,/2 ^i^p/^, 
(7a(0,0,2t,2y) (87rt/a)-'^/2e-l!^IV{2i^) ^ 

is not bounded in y G M*^. The next inequality (24) between four different 
instances of the Gaussian kernel substitutes for 3G, and so it is coined 4G. 
We note that [26, the proof of (4.4)] yields (24), too, although with a rough 
constant M (see also the first equality on p. 465 in [27] and the last equality 
on p. 15 in Friedman [12]). We also acknowledge a similar result (with 
rough constants) for the heat kernel of smooth bounded domains by Riahi 
[23, Lemma 3.1]. The optimality of the right-hand side of (24) is important in 
view of (16) and (17), and may be of independent interest. In fact, inspection 
of our calculations reveals that 6 — a in gh-a of (24) cannot be replaced by 
a bigger constant. 

Theorem 2 (4G). Let < a < b < 00 and M = f exp 

For all s < u <t and x, z,y G M"^, we have 

gb{s, X, u, z)ga{u, z, t, y) < M[gb-ais, x, u, z) V ga{u, z, t, y)] ga{s, x, t, y), (24) 
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which fails for some s < u < t, x, z,y G M"', if M < ( ) exp f -^^(^r^ 



// 1/(1 + e-i/2)< a/6 <1, then 
Proof. We have 



d/2 



d/2 



exp 



1 / \ ^1/ \ d ^ a\y — xP 
111 s, X, t, u = m47r m t — s H — ma — . 



(25) 



Considering \/2dx, v2dy and y2dz instead of x, y and z in (25), we see that 
(24) is equivalent to 

h\z — a\y — 

— ln(n — s) + m m(t — n) + m a 



< -InM + 
d 



ln(t — n) + In a 



u — s 
ln(n — s) + ln(6 — a 
a\y - 



t — u 

{b — a)\z — 



M — s 



V 



ln(t — s) + Ina 



a\y - A 

t - s 



t — u 

We rewrite this using the identity a + h — a\/ b = a f\h, and we obtain 
aU — xP r , , X , /, N (6 — a)|z — xp 



In 6 — ln(6 — o) 



+ 



M — S 



hi(ti — s) + ln(6 — a) 



u — s 



A 



ln(t — n) + In a 



a\y - A 
t — u 



2 o\v — xP 

< -lnM-ln(i- s) +lna -. 

d t — s 



Adding ln(t — u) — Ina to both sides, we have 



t - s 



+ In 



t - s 



t — u d 
a\y - z|2 



< - In M + In ■ 



b — a 



t — u 



V 



, u — s , a (5 — a)|z — xp"! 
In h In h 



t — u 



b — a 



u — s 



+ 



a\z — X 



u — s 



We denote a = a/{b — a), t = {u — s)/{t — u), E, = \y — z\\J ajit — u) and 
= \z — x\^Ja/{t — u), and observe that (t — s)/{t — u) = 1 + r. Since 
|y — x| < |z — x| + |y — z|, we see that (24) is equivalent to the following 
inequality (to hold for all r > and ^, r/ > 0), 



^^+ln(l + r)<^M + ln^+e^V 
1 + T a b 



ar 



+ In(aT) 



(26) 



We may now use Lemma 5. In fact, the constant M in (24) is optimal if 

^InM + ln^-^ = L{a). 

Considering a = a/ (J) — a) > e^^^, the last statement of the theorem follows 
from a comment after (19). □ 
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Remark 6. In applications we usually choose a (smaller than but) close 
to b, so to not loose much of Gaussian asymptotics, and in this case the 
optimality of mere M = (1 — a/h)~'^ comes as a nice feature of 4G. 



4 Discussion and examples 

In this section we discuss applications of Theorem 1 to fundamental solutions 
of second order parabolic differential operators. Namely, Theorem 1, aided 
by Theorem 2, allows for rather singular Schrodinger perturbations of such 
operators without dramatically changing the magnitude of their fundamental 
solutions. Most of the estimates given below are known, but our proofs are 
more synthetic and considerably shorter, and we have explicit constants in 
the estimates, which may be useful in homogenization problems. We also 
note that in the case of signed perturbations (not considered here) very 
precise lower bounds are obtained from Jensen's inequality for bridges [3], 
see also Remark 1 above. We begin with a discussion of Kato-type conditions 
(historical comments are given in Remark 8). 

Let d > 3. A Borel function [/ : M'' ^ R is of Kato class, if 

lim sup / dz = Q. (27) 



■5^0+ ^gKd J \Z-X 
|z-a;|<<5 



A typical example is U{z) = jzl"^''"'^, where < e < 2. By Aizenman and 
Simon [1, Theorem 4.5] or Chung and Zhao [8, Theorem 3.6], (27) holds if 
and only if for every c > the following condition is satisfied (see (22)), 



s+h 



lim sup / / gc{s,x,u, z)\U{z)\dzdu = Q . (28) 

s Rd. 

In fact, there is a constant ci = ci{d,c) such that for all /i > and U, 



\z—x\<\/h 

< ci sup / '^^ ■ (29) 

J \z — X\ 

\z—x\<\/h 

The lower bound of (29) is given in [1, (4.5)] and [8, Lemma 3.5]. The upper 
bound can be proved as in [5, Lemma 11], but for the reader's convenience 
we give a simple, explicit and more flexible argument. Firstly, by changing 
variables z = \fhz, x = ^/hx, s = hs, U{z) = hU{\/liz), we may and do 
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assume that h = 1. Then we consider the Newtonian kernel 

oo 

K{x) = j g{0,0,t,x)dt = |xp-'^r(d/2 - l)7r-'^/^/4, x £ 



nd 



and we let k{x) = g{0, 0, t, x)dt. Clearly, Jj^^j k{x)dx = 1 and k is radially 
decreasing, i.e. k{xi) > k{x2), if < \x2\- For x S M'^ and r > we 
denote B{x,r) = {y £ : \y — x\ < r}, and we consider ls(o,i), the 
indicator function of the unit ball. We observe the following auxiliary result. 

Lemma 7. // / > is radially decreasing and f is constant on B{0, 1), then 

/* 1^(0,1) > |S(0, 1/2)1/. 

Proof. We have / * Ib{o,i){x) = fB(x,i) f(y)dy- If \x\ < 1, then 

/ * lBio,i){x) > /(0)|i?(0, 1) n B{x, 1)1 > /(0)|i?(0, 1/2)1 = f{x)\B{0, 1/2)|, 

where |S(0, 1/2)| denotes the volume of B{0, 1/2). If |x| > 1, then 

/ * 1b(o,i)(^) > f{x)\B{0, \x\) n B{x, 1)1 > f{x)\B{0, 1/2)|. 

□ 

We let / = /c A A;(1,0, ... ,0) and k{x) = K{x)1b(o^i)- By Lemma 7, 

k<K*{Cdf + 6o) , 

where Cd = (|5(0, 1/2)| r(d/2 - l)7r^'^/V4)^^ and 6o denotes the Dirac 
measure at 0. The upper bound in (29) now follows for c = 1: 

\U\*k<\U\*K* {Cdf + do) < (1 + Cd) sup \U\*K. 

The case of the upper bound in (29) for general c > is similar. 

We now pass to parabolic Kato condition. For c > 0, /i > and V : 
M X R'^ M we denote 

s+h t 

-^/i(^)=sup / gc{s,x,u, z)\V{u, z)\dzdu+sup / gdu, z,t,y)\V{u, z)\dzdu. 

s,x J J t,y J J 

s Rd. t-hRd. 

We say that V is of the parabolic Kato class if lim;j^o+ ^hi^) — every 
c > 0. Considering V{s, x) = U{x) for s G M, x S W^, we may regard the 
parabolic Kato class as wider than the (time-independent) Kato class. We 
note that Nf^{V) is non-decreasing in h. Since gc is a probability transition 
density, N-,^^^^{V) < Nfi^{V) + N-,^{V) , hence 

N^_,{V)<N-^{V)+N-^{V){t-s)/h, h>0. (30) 
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To focus on nonnegative Schrodinger perturbations (in this connection see 
Remark 1), we consider, as before, Borel function q > on space-time. If 
< o < 6 < oo, then by 4G Theorem, there is a constant M' depending 
only on d and h/a, such that for ah s < i and x, y G W^, 

t 

gb{s,x,u, z)q{u, z)ga{u, z,t,y) dzdu < M'N^_^{q) ga{s,x,t,y) , (31) 

/ . \d/2 

where c = (6 — a) A a. Indeed, by (23) we can take M' = I ^ V ^ j M, 
where M is the constant from 4G Theorem. Therefore if q > belongs to 
the parabolic Kato class and < a < 6 < oo, then by (30) and (31) we have 

q(^^f{9b,ga,ib/af/^rJ,Q), (32) 

with Q{s,t) = /3(t — s) and /3 = ij/h, provided h > and > are such 
that N^l' "^^'^(q) < rj/M'. Clearly, the condition allows for arbitrarily small 
r/ > 0. 

Remark 8. The (time-independent) Kato class was first used to perturb 
the Laplace operator by Aizema and Simon [1], and was characterized as 
smallness with respect to the Laplacian on L^(M"'). The parabolic Kato class 
was proposed for the Gaussian kernel by Zhao in [26] . It was then generalized 
and used by Liskevich and Semenov [19], Liskevich, Vogt and Voigt [20] and 
Gulisashvili and van Casteren [13]. The condition is related to Miyadera 
perturbations of the semigroup of the Laplacian on L^(M'^), see Schnaubelt 
and Voigt [24]. Time- independent Kato class is wider than L'P{'R.'^) if p > 
d/2 [1], [8, Chapter 3., Example 2]. Nevertheless, the latter space is quite 
natural for perturbing Gaussian kernels, see Aronson [2], Dziubahski and 
Zienkiewicz [11], [27, Remark 1.1(b)]. Another Kato- type condition was 
introduced by Zhao in [27] to obtain strict comparability of g and g. As 
noted in [27, Remark 1.1(c)], the condition may be formulated in terms of 
Brownian bridges. This point of view was later developed in [3] (under the 
name of the relative Kato condition) and elaborated in [6] to 

t 

f f P{s,x,u,z)p{u,z,t,y) ^ J J ^ r , /QQ^ 
J J ^ ^ q{u,z)dzdu<[r] + Q{s,t)], (33) 

s X 

where s < t, x,y £ X, t] < oo, and < Q{s,u) + Q{u,t) < Q{s,t), cf. (1). 
Condition (33) indicates why we mention bridges here (see [3] for details). 
The Kato condition for bridges gives better upper bounds and seems more 
intrinsic to Schrodinger perturbations than the parabolic Kato condition, 
but the former may be cumbersome to verify in concrete situations. For the 
classical Gaussian kernel, (33) is stronger than the corresponding parabolic 
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Kato condition with a fixed c (see [3, Lemma 9] for a more general result), 
and it is rather difficult to explicitly characterize [27, Remark 1.2(a,b)]. This 
is due to the relatively large values of the integrand in (33) for (ii, z) on the 
interval connecting (s, x) and {t, y). If p satisfies 3G inequality, p{s, x, t, y) = 
p{s, y, t, x) and p is a probability transition density, then the parabolic Kato 
class and the Kato class for bridges are equivalent. This is the case for 
the transition density of the fractional Laplacian A"/^ with a € (0, 2) [3, 
Corollary 11], and the proof of this fact is similar to our application of 
4G in (31). We emphasize that each transition density p determines its 
specific Kato classes (either parabolic and for bridges), and detailed analysis 
is required to manage particularly singular q. 

Let d > 3, zi £ M^, \zi\ = 1, and B{nzi,l/n) C M'^ be the bah with 
radius \/n and center n^i, n = 1, 2, . . .. We define 

oo 
n=1 

If 5 > and n > 1/(5, then 
sup I I ^''fw-? > I n\z - nzil-'^^^ dz = f [zl^"^^^ dz . 



cglgd J \Z X\' 

\z-x\<5 B(nz-i,l/n) B{0,1) 

Therefore eU does not belong to the parabolic Kato class for any e > 0. On 
the other hand, by (29) we have, 

Nl{U) < 2ci sup / ^^%_2 dz<oo, 
xeR''- J \z — x\ 

|2-X|2<1 

hence eU G M{gb,ga, (b/O')'^^'^ ^ViQ) with linear Q and t] < 1/2, provided e 
is sufficiently small, cf. (31), (30). Since our constants are quite explicit, so 
are the resulting upper bounds for g^. To close the discussion of examples 
of q manageable by our methods, we let q{u,z) := U{z) + f{z,u), where 
/ > 0, ?7 is as above, and sup^g]gd /(u, z) is finite and increases to infinity as 
n — 7- oo. Such q requires 77 > to control U, and a superlinear Q to majorize 
sup^gffid f{u,z)du. 

We next focus on Gaussian estimates and consider a measurable transi- 
tion density p on X = 'W^ with the Borel sets Jv[ = B^d and the Lebesgue 
measure m{dz) = dz. Assume that 6 > 0, A > 1 and A G M exist such that 
for s < t we have 

p{s, X, t, y) < Ae^(*"^)gfe(s, x, t, y), x, y G M"^ . (34) 

Let <a <b and p*(s,x,t,y) = e^^^~'^^ga{s,x,t,y). In view of (23) we put 
C = A(6/a)'^/^. U q e N{gb,ga, ^, ^, which is a condition on ry and Q, 
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then q € M{p,p*,C, r], Q), and by Theorem 1, for all s < t and x,y € W^, 

Q(s.t) 

e^^'-'^ga{s,x,t,y) if0<r/<l/2, (35) 



C 



pis, X, t,y) < , 
p{s,x,t,y) < (2C)i+2Q(s,i)gA(t-s)^^(^^^^^^^)^ 



i+- 



if ?7 = . 



(36) 



In particular, such tj > may be chosen arbitrarily small, with a linear Q, if 
q is in the Kato class or the parabolic Kato class (see the discussion above) . 
This is our general setup for Gaussian estimates. 

Example 1. If p = and p* = ga, then we take A = 1, C = (b/a)'^/'^ and 
A = in (34) and, consequently, in (35,36). For clarity, a, the coefficient in 
the exponent of the Gaussian majorant, may be arbitrarily close to b, at the 
expense of the factor before the majorant in (35,36), and we require that 
q G J^igt, ga,C,rj,Q) with rj £ [0,1/2), which is satisfied, e.g., for q in the 
the Kato class. We thus recover the best results known in this setting [26, 
proof of Theorem A] with explicit control of constants. 

To relate our results to second order differential operators, we let C^(Mx]R'^) 
denote the smooth compactly supported functions on space-time, and we 
recall that for all s G M, x G and E C^(R x 



p{s,x,u,z) 



d(l){u,z) 1 . . 



dzdu 



(s,x) 



(37) 



The identity may, for instance, be obtained from integration by parts or 
by using Fourier transform in the space variable. By a general result, [6, 
Lemma 4] , the perturbed transition density p corresponds to the Schrodinger- 
type operator + g in the same way. 



p(s,x,n,z) 



+ ■^^<P{u, z) + q{u, z) 



[u,z) 



dzdu 



{s,x) 



provided q E M^gt, ga, C, Q) with t] E [0, 1/2), as above. 

Example 2. Let p be the transition density of the one-dimensional Brownian 
motion with constant unit drift, 

p{s, x,t,y) = gi{s,x - s,t,y -t) , s < t, y E M . 

There are no constants ci,C2 such that p(s,x,t,y) < cigc2{s,x,t,y) for all 
s < t and x,y £ M (cf. Zhang [25, Remark 1.3]). On the other hand, for 
each 6 E (0, 1) we have 

p{s, X, t, y) < 6^i/2g4(ib) X, t, y) , s < t, X, y E M . 

This shows why we may need A 7^ in (34), see also Norris [22]. 
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Example 3. Let / : M x M*^ 9 (s, x) ^ M, and 



Lf 



f{s,x) 



dxjdx. 



be a uniformly elliptic operator, with bounded uniformly Holder continuous 
coefficients bi and Dynkin [10, Chapter 2, l.l.A and 1.1. B] for 

detailed definitions, and [12, Chapter 1] for wider perspective. Consider the 
fundamental solution p{s, x,t,y) in the sense of [10, Theorem 1.1] for the 
following differential parabolic operator 



5/ 
ds 



+ Lf. 



(38) 



By the results of [10, Chapter 2], in particular Theorem 1.1, 1.3.1 and 1.3.3, 
p satisfies our assumptions, including (34) with — oo <ti<s<t<t2< oo. 
Therefore the Schrodinger perturbation p oi p satisfies (35) or (36) in the 
(finite) time horizon [ti,t2], if g G -^{db, da, ^, j;) V ^ [0)1/2), as 
explained before Example 1 and after the proof of Theorem 1. We thus 
recover recent results of [20, Theorem 3.10 and 3.12] (see also Remark 9 
below). We now explain how p and p are related to parabolic operators. For 
all s G M, X G M.'^ and (f> G C^{R x M.'^), we have 



p{s,x,u,z) 



d(j){u, z) 
du 



+ L(l){u, z) 



dzdu 



(39) 



In fact, the identity holds if the function cj){s,x) is bounded, supported in a 
finite time interval and uniformly Holder continuous in x, and if the same 
is true for its first derivative in time and all its derivatives up to the second 
order in space. Indeed, if we let 



h{s,x) = (l){s,x) + I p{s,x,u,z) 



d(j){u, z) 
du 



+ L(p{u,z) 



dzdu , 



then by [10, Chapter 2, 1.3.3 ] we have h = 0, which verifies (39). By (39) 
and [6, Lemma 4], the perturbed transition density p corresponds to the 
Schrodinger-type operator L + g in a similar way: for all s G M, X G R'^ and 
4> G C^{R X M'^), 



p{s,x,u,z) 



d(l){u, z) 
du 



+ L(f>{u, z) + q{u, z)4>{u, z) 



dzdu 



{s,x) 



Remark 9. In this paper by a fundamental solution we mean the negative 
of an integral inverse of a given operator acting on space-time (other authors 
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also use the terms heat kernel and Green function). More specifically, our p 
and p may be considered post-inverses of the respective differential operators 
acting on C~(M x W^), cf. Example 1 above and [27], [23, p. 13]. In the 
literature on partial differential equations it is also common to consider the 
pre-inverses, which necessitates that the applications of p, p to test functions 
be sufficiently differentiable [12, Theorem 1.5.9], [10, (1.12)]. Still differently, 
if the operator L is in the divergence form, a common notion is that of the 
weak fundamental solution, related to integration by parts, see Aronson [2, 
(2), (8), Section 6 and 7], Cho, Kim and Park [7] and Liskevich and Semenov 
[19]. It is also customary to study the operator df/ds — Lf, which is related 
to (38) by time reversal s ^ -s [10, 13], [2, (7.3)]. The setting of (38) and 
(39) is most appropriate from the probabilistic point of view: it agrees with 
time precedence and notation for (measurable) transition densities of Markov 
processes, which may be conveniently considered as integral operators on 
space-time. 

Remark 10. As we already mentioned, Zhang [27] gives sharp estimates for 
perturbations of p = g. Sharp Gaussian estimates (corresponding to p* = p) 
are generally not available by our methods if parabolic Kato condition and 
4G should be used to estimate pi. Accordingly, [27] assumes an integral 
condition related to the Brownian bridge, to bound pi by p. As we explained 
in Remark 8, the Kato condition for bridges is more restrictive than the 
parabolic Kato condition (a straightforward general approach using bridges 
is given in [6]). 

We now comment on the uniqueness of p. Trivially, p is unique because it 
is given by (2), rather than implicitly. However, in the literature of the sub- 
ject a departure point for defining p is usually one of the following Duhamel's 
(perturbation /resolvent) formulas. 



p{s,x,t,y) = p{s,x,t,y) + p{s,x,u, z)q{u, z)p{u, z,t,y)dudz, (40) 



p{s,x,t,y) = p{s,x,t,y) + p{s,x,u, z)q{u, z)p{u, z,t,y)dudz. (41) 



In short, p = p + pqp or p = p + pqp, depending on whether we consider p 
and p as pre- or post-inverses, respectively, of the corresponding differential 
operators (see [6, 4] for notation related to integral kernels). Clearly, (2) 
yields (40) and (41). Conversely, iterating (40) or (41) we get (2), and 
uniqueness, provided {pq)'^p or p{qp)^ converge to zero as n — t- oo. So is 
the case with {pq)'^p under the assumptions of Theorem 1, if p is locally in 
time majorized by a constant multiple of p*, because then {pq)"'p* — )• 0. We 
refer to [3, Theorem 2] and [20, Theorem 1.16] for further discussion of the 
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perturbation formula and uniqueness. We also note that in the setting of 
bridges, there is a natural probabilistic Feynman-Kac type formula for p [3, 
Section 6], which readily yields (2) and uniqueness. 

Here is a general argument leading to (41). We consider —p and —p as 
integral operators on space-time and post-inverses of operators L and L + q, 
respectively, which in turn act on the same given set of functions. We have 
p{L(J) + qcf)) = —(j) = pL(f), hence ptp = ptp + pqpip, where ip = Lcj). If 
the range of L uniquely determines measures, then we obtain p = p + pqp 
as integral kernels (not necessarily pointwise). This is the case, e.g., in 
the context of Example 3. We finally note that some majorization of p 
is needed for uniqueness. For instance, both p = gt and p{s,x,u, z) = 
gb{s, X, u, z) + 2dbu + |zp satisfy (37). 

5 Miscellanea 

Earlier work by Jakubowski [15] and coauthors [6] in slightly different settings 
does not require continuity assumptions on Q. Namely we call Q: M x M — )• 

[0, oo) superadditive, if 

Q{s,u) +Q{u,t) <Q{s,t) , for s<u<t. (42) 

For convenience we define Q{s,t) = if s > t. We see that t i— )• Q{s,t) is 
non-decreasing, s i— )• Q{s,t) is non-increasing, and the following limit exists, 

Q-(s,t)= lim Q(s + h,t-h). (43) 

For instance, if is a Radon measure on M and Q{s,t) = fJ,{{u G M : s < 
u < t}), then Q~{s,t) = n{{u €M: s <u < t}). 

Clearly, < Q~{s,t) < Q{s,t) and Q-{s,u) + Q-{u,t) < Q^{s,t) if 
s < u < t. We have Q = Q~ . In fact, Q~{u,v) — t- Q~{s,t) as u ^ , 
t) — 7- t~, because ifO < h < u — s < k and h < t — v < k, then Q(s + k,t — k) < 
Q~{u,v) < Q{s + h,t — h). In particular, s i— )• Q~{s,t) is right-continuous 
and 1 1—)- Q~{s,t) is left-continuous. We thus obtain the following result. 

Corollary 11. Q~{s,t) is regular superadditive. 

We note that continuous superadditive functions are used in the theory 
of rough paths by Lyons [21]. There are further similarities due to the role 
of iterated integrals in time here and in [21], and many differences related to 
the fact that we require absolute integrability (but see [16]) and also inte- 
grate/average in space (see (16) and (7)). We also note that methods similar 
to ours allow to handle gradient perturbations, which will be discussed in 
a forthcoming paper (see also [18, 16]). The next result shows that if Q is 
(plain) superadditive, then the factor rj + Q{s,t) in (16) may be replaced 
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with Cr] + CQ~ {s,t), where C comes from (12), and CQ" is regular su- 
peradditive. Thus, we may assure regular superadditivity at the expense of 
increasing rj and Q. 

Lemma 12. Assume that p and p* are transition densities, function q>0 
is measurable on space-time, rj > 0, C > 1, Q is superadditive, and (12) and 
(16) hold. Then for all s < t < oo and x,y £ X , we have 

s X 

Proof. For x,y £ X , s < t, h > 0, hy Chapman-Kolmogorov, (16) and (12), 

t-h 

p{s, x,u, z)q{u, z)p* (u, z,t,y) dzdu = J J p{s, x, s + h,v) 

s+h X XX 
t-h 




p{s + h, V, u, z)q{u, z)p* (n, z,t — h, w)dzdup* {t — h, w, t, y) dwdv 

s+h X 

< \i] + Q{s + h.,t — /i)] j p{s, x,s + h, w)p* (s + h, w, t, y) dw (44) 

X 

< C[r] + Q{s + h,t - h)]p*{s,x,t,y) . 

We then let /i — )• 0, and use (43) and the monotone convergence theorem. □ 

If p* is a time-changed p, then we can do even better. 
Lemma 13. Under the assumptions of Lemma 12 we have 



t 



p{s,x,u, z)q{u, z)p* {u, z,t,y) dz du < [r] + Q {s,t)]p* {s, x,t,y) . 



s X 



if{s, oo) 9 1 1— 7- p(s, X, t, y) is continuous, p is time-homogeneous: p{s, x, t, y) = 
p{s + r,x,t-\-r,y) forr £ W, and p*{s,x,t,y) = p{as,x,at,y) for some a > 0. 

Proof. Picking up (44), for s < t, x,y ^W^, we have 

lim sup / p{s, x,t — h, w)p* {t — h, w, t, y) dw 
/1-S.0+ J 

= lim sup / p{s, x,t — h, w)p{t — h,w,t ah — h, y) dw 

h~^0+ J 
Rd 

= limsupp(s, x,t -\- ah — h,y) = p{s, x, t, y) . 

/l-5>0+ 

□ 
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Lemma 13 applies to the Gaussian density, if p = gj, and p* = ga, where 
< a < 6. Indeed, gf,{s,x,t,y) = ga{as/b,x,at/b,y), see also (23). 
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